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Abstract 

We investigate to what extent the experimental information on i? — > vrvr branching 
fractions and CP asymmetries can be used to better understand the QCD dynamics in 
these decays. For this purpose we decompose the independent isospin amphtudes into 
factorizable and non-factorizable contributions. The former can be estimated within 
the framework of QCD factorization for exclusive B decays. The latter vanish in the 
heavy-quark limit, mf, — > oo, and are treated as unknown hadronic parameters. We 
discuss at some length in which way the non-factorizable contributions are treated 
in different theoretical and phenomenological frameworks. We point out the potential 
differences between the phenomenological treatment of power-corrections in the "BBNS 
approach" , and the appearance of power-suppressed operators in soft-collinear effective 
theory (SCET). On that basis we define a handful of different (but generic) scenarios 
where the non-factorizable part of isospin amplitudes is parametrized in terms of three or 
four unknowns, which can be constrained by data. We also give some short discussion 
on the implications of our analysis for B — > ttK decays. In particular, since non- 
factorizable QCD effects in B ^ vrvr may be large, we cannot exclude sizeable non- 
factorizable effects, which violate SU{3)f flavour symmetry, or even isospin symmetry 
(via long-distance QED effects). This may help to explain certain puzzles in connection 
with isospin-violating observables in B ^ nK decays. 



^Heisenberg Fellow 



1 Introduction 



An important unresolved question in the theoretical analysis of charmless non-leptonic B 
decays is the quantitative understanding of the non-perturbative dynamics responsible for 
non-factorizable contributions to decay amplitudes. 

In the heavy quark limit, rrif, — > oo , a factorization theorem pP states that B PP 
decay amplitudes (where P stands for a light pseudoscalar meson) factorize into pertur- 
batively calculable coefficient functions Tjji and universal hadronic quantities, namely 
B —> P transition form factors F^^^{mj>) and light-cone wave functions for heavy and 
light mesons (0^ and (f).^). Schematically, one has 

{PP\H^s\B) = F^^^ ■Ti0(Pp + Tu ®<Pb®<Pp<^<Pp 

+ terms suppressed by l/m^ (1) 

where the symbol ® represents convolution with respect to the light-cone momentum frac- 
tions of light quarks inside the mesons. Without taking into account radiative corrections, 
Ti depends on kinematic factors only, and Tn vanishes. This approximation corresponds 
to the "naive" factorization assumption. At first order of the strong coupling constant 
the factorization formula has been shown to hold by the explicit calculation of the 
corrections to naive factorization [H I2] • In the following we refer to in the heavy-quark 
limit as "QCD factorization" . Arguments towards an all-order proof have been given in an 
effective theory framework in [Hj (see also ^ for a recent discussion). Here, the relevant 
momentum modes for B — > PP decays are determined by the momentum scaling of the 
external particles and their possible interactions. The QCD factorization formula follows 
from identifying Tjji as coefficient functions of operators in the soft-coUinear effective the- 
ory (SCET), where the constituents from different hadrons appear to be decoupled. An 
all-order proof, which considers the subtle effects of endpoint divergences (see Section 2 in 
[5] for a toy example, and also |5]), has not been worked out to the last detail. It should 
follow a similar line of reasoning as for the somewhat simpler cases of -B — >^ 7 and B ^ ir 
form factors [3 01312] • 

A major complication for phenomenology arises due to the fact that at least some of 
the power corrections to appear to be enhanced by large numerical coefficients (these 
are proportional to the quark condensate in QCD, for this reason these terms are referred 
to as "chirally enhanced"). In the diagrammatic analysis, non-factorizable contributions 
are identified from convolution integrals that suffer from endpoint- divergences when one 
of the parton energies vanishes. The authors of [1, parametrize the chirally enhanced 
contributions of this type by an arbitrary complex number. The modulus of this num- 
ber is estimated from regularizing the endpoint divergences with a finite energy cut-off. 
We emphasize that this is a mo del- dependent procedure that aims to get a quantitative 
handle on terms that are beyond the QCD factorization approach. We will refer to this 
approximation and its phenomenological implications as the "BBNS approach" . 

Another well-known framework is to use approximate flavour symmetries (isospin or 
SU{3)f) to relate different decay amplitudes and reduce the number of unknown hadronic 
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parameters P CD] • This procedure is often combined with so-called "plausible dynamical 
assumptions" about the importance of certain flavour topologies that can be identified in 
the factorization approximation only. In particular, from a recent study in ^l] along these 
lines, it has been concluded that (a) non-factorizable effects in 5 — »• tttt are large, (b) 
certain B nK observables point to possible new physics effects in electroweak penguin 
contributions. 

A comprehensive analysis of presently available B PP data leads to a somewhat 
milder conclusion It is found that, within the statistical uncertainties, the "anomalies" 
in i? — i> ttK decays may still be considered as being consistent with the SM. Furthermore, 
the non-factorizable effects in i? — vrvr can be more or less accounted for by fitting the 
hadronic input parameters in the BBNS approach to experimental data. Of course, in this 
procedure, most of the predictive power of the QCD factorization approach is lost. Still, 
the model-dependent parametrization of non-factorizable effects in the BBNS framework 
turns out to result in significant constraints when used as the basis for a CKM fit. 

The purpose of this article is to carefully examine different (model-dependent) ap- 
proaches to quantify non-factorizable hadronic effects in charmless non-leptonic B de- 
cays. In view of the ultimate goal, namely to extract independent information on CKM 
parameters from non-leptonic B decays, we think that it is important to make sure that 
certain assumptions about the size of flavour symmetry breaking, the origin of strong 
phases etc. are clearly identified in order not to induce an uncontrolled theoretical bias. 

Our paper is organized as follows. In the next section we will use the 5 — tttt de- 
cays as a guideline, and decompose the independent isospin amplitudes into factorizable 
and non-factorizable parts. The factorizable contributions are estimated within the QCD 
factorization approach to first order in the strong coupling constant, and using default val- 
ues for hadronic input parameters. The non-factorizable contributions are considered as 
unknown parameters, the size of which has to be taken from experimental data on branch- 
ing fractions and CP asymmetries. We will discuss the origin/interpretation of different 
sources for non-factorizable effects within the BBNS approach, SCET, and phenomeno- 
logical studies assuming large long-distance penguin contributions. From this we develop 
a handful of constrained scenarios that implement different generic features of such ap- 
proaches. Constraining the parameters for these scenarios using B ^ data, one may 
obtain quite different results on the size of individual isospin amplitudes. In particular, 
we find that in scenarios with large non-factorizable penguin or annihilation contributions, 
one may not exclude sizeable corrections to isospin-violating observables which arise from 
long-distance QED effects, and may be relevant to explain the B — > nK puzzles men- 
tioned above. A detailed numerical analysis of i? ^ nK modes is postponed until new 
experimental data become available. 
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2 Factorizable and non-factorizable contributions 



2.1 Isospin decomposition of decay amplitudes 

The CKM elements entering the effective weak hamihonian for 6-quark decays in the 
Standard Model are defined as 

aS'^ = v^im:, ■ (2) 

In particular, for b d transitions we have 

Ac 

and for b —>■ s transitions we have 

\(«) 

tan'' Oc Ruc-'^ = eKMe-'^ (4) 



A, 



is) 



with sin^c" — 0.2266 being the Cabbibo angle, and Ru = ^/p' + fj'^ — 0.405 and 7 ~ 
62° related to one side and one angle of the CKM-triangle (all numbers from p!^, i?„ is 
sometimes denoted as R^). In the case of 6 — >■ d decays, all CKM factors are of the same 
order, whereas for 6 — > s decays, e^M is suppressed by two powers of the Cabbibo angle. 

Assuming isospin conservation in hadronic matrix elements the B —* ittt decay ampli- 
tudes can be decomposed into^ 

V2{7r-7r'\H,s\B-) ^ X^f^ [3A„(2,3/2)] + A^^^) [3A,(2,3/2)] , (5) 

(7r+7r-|i7eff|5°) ~ AW[-A„(0,1/2) + A,(2,3/2)] , 

+A('^)[-A,(0,l/2)+A,(2,3/2)] (6) 

v^(7rV|i7eff|5°) ^ AW[A„(0,1/2) + 2A„(2,3/2)] 

+A('^)[A,(0,l/2) + 2A,(2,3/2)] , (7) 

where the first argument denotes the total isospin / of the final state, and the second 
argument denotes the isospin A/ of the operators in the weak effective hamiltonian. For 
the charge-conjugated modes one has to replace A-'^'' by A-'^''*. Thus, the physical B tttt 
amplitudes fulfill the well-known isospin relation 

A[B- ^ 7r-7r°] - ^ 7r°7r°] - ^ tt+tt-]/V2 = . (8) 

The latter is violated by small quark mass effects and electromagnetic corrections which 
we discard in the following. 



^The |7r°7r°) state in this notation already includes a statistical factor 1/V2 from Bosc symmetry, i.e. the 
branching ratio calculated with this amplitude does not receive an additional factor 1/2. This corresponds 
to the convention in and differs from the convention in 
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Notice that only electroweak penguin operators contribute to the amphtude Ac (2, 3/2). 
These operators have very small Wilson coefficients compared to those in ^4^(2, 3/2). Ac- 
tually, neglecting the tiny Wilson coefficients C^ and Cs, one can use Fierz identities to 
relate the matrix elements of ^9,10 with those of Oi^2 and thus obtains 

M'^.'i/'^) _ 3 + Cio _ 

A„(2,3/2) - 2 + ^ ■ ^ ^ 

The contribution from Ac(2,3/2) can therefore be neglected in i? — tttt decays to a very 
good approximation. This result is based on the structure of the effective electroweak 
hamiltonian and on the general isospin analysis only ITH] . 

Using this simplification one sometimes introduces an intuitive parametrization, which 
refers to the flavour topology ("penguin", "tree" or "colour-suppressed") of amplitudes 
that one would obtain in the factorization approximation, 

P = A('^)*Ae(0,l/2) , 
Te^^ = Af*(A„(0,l/2)-A„(2,3/2)) , 

e'en = -AW*(A,(0,1/2) + 2A„(2,3/2)) . (10) 

The ratios between these amplitudes can be parametrized in terms of two moduli and two 
strong phases as 

de-^-^e- = -^-^Q'V^) (11) 

f A„(2,3/2)-A„(0,l/2) ^ ' 

as a measure for the "penguin-to-tree ratio" , and 

^ £ ^ A„(0,l/2) + 2A.(2,3/2) 

-f A„(2,3/2)-A.(0,l/2) ^ ^ 

as a measure for the ratio of "colour-suppressed" to "colour-allowed" tree amplitudes. 
Beyond the factorization approximation the notion of flavour topologies might be somewhat 
misleading, whereas the classification in terms of isospin amplitudes is more general. 

The isospin amplitudes Aj(/, A/) contain contributions from short-distance dynamics 
(modes with large virtualities in the heavy quark limit) and long-distance dynamics (modes 
with virtualities of order Aq(jj-)). The short-distance effects can be treated in perturba- 
tive QCD, making use of the heavy-quark expansion. The long-distance effects represent 
hadronic uncertainties. We therefore decompose every isospin amplitude as 

A,(J, AJ) = (/, A/) + (/, A/) . (13) 

In "naive" factorization, the amplitudes Af can be expressed in terms of electroweak 
Wilson coefficients and hadronic decay constants and form factors. An example for a 
naively factorizing diagram is shown in Fig. where we also indicate the momentum 
scaling of external and internal lines Here and in the following, "s" stands for soft 
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Figure 1: Example for a naively factorizing contribution to B —>■ PP. The labels indicate 
the scaling of the momentum modes. 

momenta, "c" and "c" stand for coUinear momenta (virtuality A^) in one or the other 
direction, "he" and he" for hard-collinear momenta (virtuality Am^,) in one or the other 
direction, and "h" for hard modes. 

In the heavy-quark limit, we can use QCD factorization (Q) to improve the quantitative 
description of the factorizable part A[ . The contributions to the first term in (Q) come from 
vertex corrections and penguin contractions as shown in Fig. and from hard spectator 
scattering as shown in Fig. El^a). Non-factorizable effects arise from power corrections in 
the Aqcd/'"^6 expansion. They include the annihilation topologies (see Fig. H)), and cannot 
be calculated in a reliable way at present. 



Figure 2: Examples for (a) vertex and (b) penguin corrections to naive factorization in 
B PP. The labels indicate the scaling of the momentum modes. 



2.2 Factorizable contributions 

In |rj the contributions related to Tj and Tn in are expressed in terms of parameters 
Qi = cii,i + CLi,u (yCLi,!! will be restricted to the heavy-quark limit, whereas in ajj the terms 
proportional to = m^/mqmiy are kept). The factorizable amplitudes Af read 
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(14) 



which determines the overall normalization of the factorizable amplitudes, and F( 
is the scalar B n transition form factor. 




Figure 3: Hard-scattering contributions to i? ^ PP. (a) Example for one-gluon exchange 
included in the BBNS analysis, (b) Example for (power-suppressed) higher-order diagram. 
The labels indicate the scaling of the momentum modes. 

The error in the first term of ()14j) refers to the variation of the factorization scale fi 
between mb/2 and 2mb in the vertex and penguin graphs. The second term denotes the 
central value for the hard-scattering contribution which has a large uncertainty related 
to the first inverse moment A^^ of the light-cone distribution amplitude of the B meson. 
(There are more sources of parametric uncertainties, in particular the scale-dependence of 
the hard-scattering term, see the numerical discussion in |lj|. Notice that for the hard- 
scattering terms, we considered the electroweak Wilson coefficients at the scale /i = nib-) 
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(b) 
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Figure 4: Annihilation contributions to -B — > PP. (a) Example for one-gluon exchange 
included in the BBNS analysis, (b) Example for higher Fock-state-contribution. 

In the following numerical discussion we will refer to the central values of the A[ only. 
Of course, a significant change in the numerical values of Af would also influence our con- 
clusions about the size of Af^^. In particular, one has to keep in mind that quantities like 
the "colour-suppressed" amplitude C defined in fllO|) involve large cancellations in naive 
factorization, and therefore are particularly sensitive to the precise value of hard-scattering 
and non-factorizable corrections. The question whether the variation of all possible input 
parameters in the BBNS approach within "reasonable" ranges could reproduce the exper- 
imental data has already been studied in |T2I.^ Here we take the point of view that the 
default values in ^ give a reliable prediction for the factorizable contributions, whereas 
the BBNS analysis of non-factorizable terms (including the variation of hadronic input 
parameters) is considered as one option among different alternatives. 

The power-corrections to the hard-scattering parameters ajji, and the annihilation pa- 
rameters bi are considered as part of the unknown functions A^^, which we parametrize 
as 



with r, > and an arbitrary phase 0,. (As explained above, one can safely set A^^{2, 3/2) = 
0: Even if we allow for an order-of-magnitude enhancement with respect to its factorizable 
counterpart in we would only get an 0{1%) correction to ^4^(2, 3/2).) 

2.3 Non-factorizable effects from B ^ tttt data 

Our general parametrization of non-factorizable effects introduces seven adjustable para- 



meters r,(0,l/2), 0^0,1/2), r„(2,3/2), 0„(2,3/2), r,(0,l/2), 0e(O,l/2), and Fo^^\ On 



the other hand, if we neglect the tiny contribution from Ac{2,3/2), on general grounds, 
we have only five relevant parameters to describe three complex isospin amplitudes for 
B — »• TTTT (one overall phase is not observable). Consequently, our parametrization contains 
some redundancy, which we will keep for the moment. Later we will consider different 
constrained scenarios, where the number of parameters is less than 5. 

^For instance, using the "large-a2" scenario in 0, where (u^^)-^ = 4.2, As = 200 MeV, and F^^^'" = 
0.25, the hard-scattering correction in aiji increase by a factor of 3. 



Af^(J,AJ) := n(J,AJ)e*^'(^'^^) A 



(15) 
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Observable 


BaBar 


Belle 


CLEO 


Average 


B[B^ 71 + 71-] 

B[B+ tt+ttO] 


4.7 ±0.6 ±0.2 

5.5tJ:|]±0.6 
2.1 ±0.6 ±0.3 


4.4 ±0.6 ±0.3 
5.0 ±1.2 ±0.5 
1.7 ±0.6 ±0.2 


A e-+1.4+0.5 
^•^-1.2-0.4 

^•0-1.6-0.7 

< 4.4 


4.55 ±0.44 
5 1 8+°-'^'^ 

^■-'-°-0.76 

1.90 ±0.47 




-0.19 ±0.19 ±0.05 
-0.40 ±0.22 ±0.03 


-0.58 ±0.15 ±0.07 
-1.00 ±0.21 ±0.07 




-0.46 ±0.13 
-0.73 ±0.16 


^cp[vr+7r0] 


-0.03l^;{? ± 0.02 


-0.14 ±0.24^^:^^ 




-0.07 ±0.14 



Table 1: Experimental results on i? ^ tttt observables [12]. All branching ratios are CP- 
averaged and quoted in units of 10"^. The sign convention for CP-asymmetries is defined 
in the text. 

Our strategy to infer information on the non-factorizable parameters from experimen- 
tal data is to produce sets of random parameter values, and calculate the corresponding 
X^-value by comparing the theoretical branching ratios and CP asymmetries with the ex- 
perimental measurements in Table We follow the sign convention of |12j . 

l[B ^ TT TT^J ± 1 [B^ 7l + 7l^\ 

and 

C+- ^ 2ImA^7r ^ 1 ~ IAttttP 

Ti-TT 1i|\|2' 1_|_|X|2' ^ ' 

q A[B^^ 7^+7.-] ^ _,., ^[E°-.vr+vr-] 

The so-obtained distributions enable us to investigate the generic size and importance 
of non-factorizable parameters. To generate the sample, we assume uniform distributions 
of parameter values in the following ranges 

0.23 < F^^^{ml) < 0.33 , 
0<r,,,(/,A/) < 1.0 , 
0° < 0„,e(/, A/) < 360° , (19) 

The bound on the scalar form factor is the main theoretical bias. We used a rather 
conservative estimate of the theoretical uncertainties, which contains the central value 
0.28 used in [1. as well as a recent update 0.26 for this quantity in the framework of QCD 

''Updated results for charmless non-leptonic B-decays have been presented by the BaBar and Belle 
collaborations in jlHI 1141 IT^ llfij . The new data are consistent with the numbers presented in jT2]. We 
have checked that they would only lead to marginal effects in our numerical analysis, and our qualitative 
conclusions remain unchanged. With the new data the compatibility with (|2DJ), see below, appears to be 
slightly improved. 
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sum rules [12]. The upper bound on ru,c{I, A/) will turn out to be sufficiently large not to 
induce an additional bias. 

Already for the unconstrained scenario we find some interesting patterns, see Fig. 
where we have plotted the value against each of the free parameters for a sample of 500 
points with < 10: 

• The value of the scalar B ^ it form factor is not constrained by the data^ 

• The distribution of the parameter r„(0, 1/2) is rather broad, with a slight prefer- 
ence for r„(0,l/2) = 0.5 ±0.5. The corresponding phase take values in the range 
120° < (j)u{0, 1/2) < 240°. (Notice that for values of 0^(0, 1/2) around 180° the non- 
factorizable contributions reduce the value of the "colour-allowed" tree amplitude, 
which accomodates the experimental fact that the parameter x in ()12|) is large.) 

• We further find rc(0,l/2) < 0.35, and r„(2,3/2) < 0.8. The corresponding phases 
0c(0, 1/2) and 0m(2, 1/2) show no pronounced preference. 

Of course, there are correlations between the different parameters. For instance, there 
is no solution with small x^? where all parameters ri{I,AI) are small (relative to the fac- 
torizable terms), in other words the default scenario in the BBNS approach is disfavoured 
by the data (in accordance with similar conclusions in |1H IT^ 1^). For illustration we 
marked solutions which simultaneously fulfill 

r„(0, 1/2) < 0.5 , r„(2, 3/2) < 0.2 , and r^{0, 1/2) < 0.1 , (20) 

and which, in view of the large parametric uncertainties related to hard-scattering and 
annihilation contributions (see below) may still be viewed as more or less compatible with 
the BBNS approach (in the spirit of ll2j). This still yields values of — 3. 

2.4 Non-factorizable effects from "BBNS" 

In the BBNS approach non-factorizable effects arise through chirally enhanced power- 
corrections which are identified from endpoint-divergent convolution integrals, appearing 
in the diagrammatic approach. 

2.4.1 Hard-scattering contributions 

One source of non-factorizable power-corrections are so-called hard spectator-scattering 
diagrams, where a gluon connects the spectator quark to the short- distance decay process. 
Fig. ISfa) (together with an analogous diagram where the gluon is attached to the other 
collinear quark) has been considered in ^ . On the other hand. Fig. Efb) represents an 
example of a higher-order diagram (which is not included in the BBNS approach), which 

^It would therefore be interesting to independently measure the value of the form factor from B — > tt£i/ 
decay for given value of \Vub\- 
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Figure 5: Comparison of values for a random sample of non-factorizable parameters 
Dots for unconstrained scenario, stars for combinations that fulfill ((201) • 
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involves a multi-particle Fock state in the B meson and which is power-suppressed (see 
below) . 

Apart from a factorizable part that determines the heavy quark limit, explicit calcula- 
tion shows that the diagram in Fig. Efa) gives rise to chirally-enhanced power-suppressed 
endpoint divergences. At the considered order in the diagrammatic expansion, these end- 
point divergences enter through the quantity 



l-Afe/ms 



U 



U 



■4>l{u) 



(21) 



where (ppiu) ~ 1 is a twist-3 light-cone distribution amplitude of the pion. We regularized 
the integral by means of an ad-hoc cut-off, because the integral does not converge for 
M — s> 1. Since the twist-3 distribution amplitudes are normalized to a term proportional to 
the quark condensate (i.e. the ratio of Goldstone boson and quark masses, m\^/mg), their 
contributions are numerically large ( "chirally enhanced" ) . To study the phenomenological 
impact of these terms, the authors of Ij propose to parametrize the quantity on the left- 
hand side in terms of 



'1 + ph e 



In 



rriB 



A 



H 



PH ~ 0(1] 



(22) 



with ifH being an arbitrary phase, whereas Ah ~ 0.5 GeV. 

Inserting the default values for the hadronic parameters in the BBNS analysis, we 
obtain 



^r[2,3/2]ks 



-OmOA^^{l+ Ph e'^» 
-0.018 A,^(1+Pj^ e'^P" 
A^^'[0,l/2]|h. ~ 0.001 A^^(l + pHe''^«) 



0.28 350 MeV (m^^)^' 




\b 3.3 
0.28 350 MeV {u-^)^ 



^0 



A 



B 



3.3 



Ai^^[2,3/2]|h. ^ 0. 



(23) 

(24) 

(25) 
(26) 



Comparing with the general parametrization (fT3j) . we deduce that in the BBNS approach 
r„(/. A/) are expected to receive contributions of the order 5% (up to 10% for the "large 02 
scenario") from hard spectator-scattering, whereas the corresponding effect for rc(/. A/) 
seems to be negligible. 

In any case the above procedure is understood to only give a rough idea about the 
typical size of non-factorizable contributions for individual decay amplitudes. The origin 
and size of strong re-scattering phases remains unclear. In particular, treating pn and ipn as 
universal parameters, one induces model-dependent correlations between non-factorizable 
effects in different isospin amplitudes. Substantially different strong interaction phases for, 
say, / = and 1 = 2 final states can arise from more complicated diagrams with additional 
quark lines (i.e. higher Fock states). In the diagrammatic approach along the lines of BBNS, 
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these effects can only sliow up at liiglier orders in tlie diagrammatic expansion. Notice tliat 
in tlie case of non-factorizable endpoint configurations, "liiglier-order" diagrams are not 
necessarily suppressed by powers of the strong coupling constant. The universality of strong 
phases arising from hard-spectator scattering is therefore not a generic feature of QCD 
factorization, but represents an additional phenomenological assumption in the BBNS 
approach. The analysis ^2] uses universal BBNS parameters for all isospin amplitudes. 
Different phases in different isospin amplitudes can still be generated by combining the 
non-factorizing pieces from hard scattering (Xh) and from annihilation [X^ see below) 
which are treated as independent complex parameters. Nevertheless, a model-dependent 
correlation between non-factorizable contributions to different B ^ tttt and B ttK 
amplitudes remains. The implications of that analysis should therefore be interpreted with 
some care. 

2.4.2 Scenario 1: Dominance of hard spectator scattering 

If we assume that non-factorizable effects from power-suppressed contributions to hard 
spectator scattering (or alternatively, significant changes of hadronic input parameters in 
the factorizable pieces) are the main source for A^^{I,AI), we may define a phenomeno- 
logical approximation, where we put rc(0, 1/2) = 0. In addition, we might either fix the 
form factor value, Fq = 0.26 (Scenario la), or assume that the phases (pn are universal 
such that 0„(0, 1/2) = 0„(2, 1/2) (Scenario lb). 

Repeating the analysis of the B tttt data with these additional constraints, we obtain 
the situation illustrated in Figs. IHl and [7| Scenario la still gives a more or less reasonable 
description with ^ 4 for six experimental observables and four adjustable parameters. 
The situation in Scenario lb is similar, but somewhat worse. In both cases we need 
rather large non-factorizable amplitudes with either ru(0, 1/2) > 0.5 or r^j(2,3/2) > 0.3. 
0u(0, 1/2) is rather constrained in both cases. 

2.4.3 Annihilation topologies 

Another source of non-factorizable, power-suppressed contributions to B ^ PP are an- 
nihilation topologies, see Fig. |31 They receive the same chiral enhancement as the non- 
factorizable hard-scattering pieces. 

In the BBNS approach, the flavour- changing sub-process for annihilation topologies is 
bq qq, and consequently, it can only contribute to / = (for B — » vrvr) or I = 1/2 
(for B ttK) decay amplitudes (again this statement is only true as long as we do not 
consider higher Fock states). The order of magnitude for annihilation effects is estimated 
in a similar way as for the hard-scattering terms in the previous section, introducing the 
quantity 

AI = (l + p^e*^-)ln^, PA~0(1). (27) 

The corresponding contributions to the non-factorizable B ^ irn amplitudes read 

A„(0, l/2)Unn ^ 0.01 (1 + l.Op^e^'^-* - O.Tp^e^^^-^) A^^ , (28) 
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Ae(0, l/2)|ann ^ -0.01 (1 + 2.0pAe''^^ + O.gpie^^^-^) A^^ . (29) 

Comparing with the general parametrization ^T^ . we deduce that in the BBNS approach 
r^i c[0, 1/2] may receive contributions of the order of several percent from annihilation. 
Notice that the annihilation topologies in BBNS lead to l/ml corrections that are doubly 
chirally enhanced (giving rise to the p\ terms above). 

2.4.4 Scenario 2: Dominance of annihilation topologies 

In another approximation, we assume that the non-factorizable part of the A J = 3/2 
amplitude is sub-dominant and can be neglected, r„(2,3/2) = 0. In addition we again fix 
the form factor, F^^"" = 0.26. 

From the theoretical point of view, the situation corresponds to the case where the 
annihilation graphs in the BBNS approximation are assumed to be the dominant source 
of non-factorizable effect. Alternatively, it can also be viewed as representing the case 
where higher-order contributions from penguin corrections ( "charm" and "GIM" penguins 
[^122]; for a recent phenomenological fit along these lines, see |23|) are the main source 
of non-factorizing effects (again, such effects do not contribute to the amplitudes with 
AI = 3/2). 

From the plots in Fig. |Hlwe observe that a very good description of the data is possible 
in such a scenario. Again, rather large non-factorizable effects (i.e. of the same order as 
the factorizable ones) are needed in ^4^(0, 1/2) and/or Ac{0, 1/2). The data also show a 
clear prefererence for ^^^(0, 1/2) < A^^{0, 1/2) in this scenario. 

2.5 Non-factorizable contributions in SCET 

Soft-collinear effective theory has been developped as a systematic tool to study the factor- 
ization of different short- and long-distance modes contributing to inclusive and exclusive 
B decays j21l 123 I2S1 • In case of exclusive decays, two kinds of short-distance modes are 
successively integrated out: hard modes (with virtualities of order ml) and hard-collinear 
modes (with energies of order m,b and virtualities of order mbAqcD) [2Z1I2H1IHI- The first 
matching step leads from QCD to the so-called SCETj. The second matching step leads 
from SCETi to SCETn. The effective theory SCETn only contains long-distance modes 
with virtualities of the order of the QCD scale. In the B meson rest frame these are 
denoted as "soft" (all momentum components of order Aqcd)? and "collinear" (one mo- 
mentum component scales as rrib)- Fields and operators in the effective theory have a 
definite power-counting in terms of the expansion parameter = Kqcd/itii,. In this para- 
graph we will discuss some generic examples for effective-theory operators that may give 
rise to factorizable and/or non-factorizable effects in B ^ mr. 

Let us follow the strategy of [5^ (which has been used in the context of a factorization 
proof for the B Triu decay) and identify the possible field content of effective operators 
in SCETi that are relevant for B — >■ tttt. (We will consider light-cone gauge for the collinear 
modes, and will drop soft Wilson lines for simplicity. We also do not explicitely note Dirac 
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or Lorentz indices.). The two hard-coUinear directions defined by the final-state hadron are 
denoted as "he" and "he" , respectively. It is understood that for all operators that we will 
list below, one has a corresponding term with he ^ he interchanged. It should also be clear 
that the different possible flavour structures of the operators in SCET are obtained from 
matching the corresponding operators in the weak effective hamiltonian by integrating out 
hard QCD modes (we will comment on flavor- dependent QED effects below). 
In the second matching step, one has to generate the minimal fleld content 

[qs... K] [^c • ■ ■ ^c] [^c ■ • • ^c] 

that is necessary to build up the initial and flnal state quantum numbers (the dots stand 
for additional qq pairs or gluon flelds of the same kind; we do not consider decays into 
flavor-singlet mesons here). The generation of soft and coUinear flelds from hard-collinear 
ones costs a certain power of the small expansion parameter A, which can be read off the 
corresponding interaction terms in the SCETj Lagrangian. Examples are j3] 



^hc 


2^ 




^hc 






^hc 






^hc 







2.5.1 3-body operators 

The minimal possible fleld content for a SCETi operator, that contributes to B tttt, is 

[ehcA^/^.]~A^ (31) 

where on the right-hand side we indicated the power-counting for this operator, following 
from the SCETj Lagrangian. Performing the explicit matching calculation, one flnds that 
at tree-level the flrst non-trivial operator is just a copy of the chromomagnetic term Of in 
the weak hamiltonian, restricted to the particular kinematical situation,^ 

-Cf 1^ 4c a^.(l + 75) K + (he ^ h^) . (32) 

Here Cf^ contains the short- distance contribution from loop contractions with 4-quark 
penguin operators. 

Starting from the operator in (jH^ . in order to generate the necessary flnal-state partons, 
we need at least two coUinear quark flelds C,c and C,c which have to be generated from the 
hard-collinear fleld C,hc, and two coUinear quark flelds and which should descend from 
A^. The flrst case costs at least a factor A'^, for instance through the chain 

^cAt ^ UUc] . (33) 

^The power-counting follows from the leading term n'LG^^a = {Ti^d)A^ ^ A in the light-cone gauge. 
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In the second case, we cannot directly produce two coUinear quark fields from the hard- 
collinear gluon field, because they would be in a flavor- singlet configuration (the case of 
flavor- singlet mesons has been discussed in the context of QCD factorization in jSH]). There- 
fore we need at least two additional quark fields that do not end up in the corresponding 
pion (and therefore have to come from the initial B meson which provides soft modes). A 
possible branching 

^Ec^^.^h^^fcediesgs] (34) 

costs a factor such that the power-counting for currents in SCETn that descend from 
three-body operators in SCETi is A^"^. Together with the power-counting for the hadronic 
states the contribution to the B tttt amplitude is A^, which has to be compared with 
the result in naive factorization ~ J^^Fq^'^ ~ A^. Therefore, contributions from three- 
body operators are l/mt suppressed (and therefore do not contribute to the factorization 
theorem in the heavy-quark limit). To obtain a non- vanishing contribution of a 3-body 
SCETi operator to a 5 — tttt matrix element in the diagrammatic expansion of the 
BBNS approach, we thus need at least three gluons, see Fig. Ofb) . Such diagrams are not 
considered in the analysis of power-corrections in PP, which has been restricted to one-gluon 
exchange diagrams. Up to now, we do not know whether this operator gives factorizable 
contributions in the first non- vanishing order in SCETn (i.e. 1/mf, with respect to the 
leading contributions from naive factorization). In any case, from the structure of the 
factorization proof for B —>■ iriu decays in and similar arguments in [32], we expect 
that, in general, factorization of soft and coUinear modes in SCETn does not hold for 
power-corrections obtained in the matching of SCETj to SCETn. Therefore, at least on 
the level of 1/ml power-corrections, the operator in (j32|) provides a new source of non- 
factorizable corrections. Similarly as for the annihilation diagrams considered in BBNS, 
see (j2Z|), contributions from such power-suppressed 8-quark operators in SCETn could 
be doubly chirally enhanced. In this case, numerically they may be as important as the 
non-factorizable terms included in the BBNS analysis. 

2.5.2 4-body operators 

The leading-power contributions to the function Tj in (JJ) involve four-quark operators in 
SCETi of the type 

l^T^^i^MhcK] \' . (35) 

The conversion of the hard-collinear quark pair to a coUinear one costs twice a power of 
A, and the conversion of via (j3^ costs a factor of A"^. Therefore, these four-quark 
operators match onto 6-quark operators in SCETn which scale as A^^ and, according to 
the discussion in the previous paragraph, are leading power. 

Another 4-quark operator that is allowed in SCETj by momentum conservation is of 
the form 

fc^.]fehc] ~A« . (36) 
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The conversion of hard-collinear quark fields into at least two collinear fields via ()33|) 
costs each a factor of A'^. Therefore, these operators match onto 8-quark operators in 
SCETii which scale at least as A^^ and are thus suppressed. This corresponds to 

an annihilation topology which involves a higher Fock state {bq){q'q') in the B meson. 
Notice that, again, the suppression by 1/ml can be numerically compensated by two chiral 
enhancement factors coming from the wave functions of the two final-state mesons. 
We may also consider 4-body operators that involve additional gluons, like 

[A^] [^^.Ai^K] ~ A^ J... im l^-cQs] [Uc] - , (37) 

or 

[AjAiM-^K]r^\' I . . . Aj[UM^Qs][qsmcK] - (38) 

where again we have to use (jMj) to obtain flavour non-singlet collinear quark configurations 
from a single hard-collinear gluon field. 

2.5.3 Remark on the treatment of charm quarks 

One may also worry about 4-quark operators involving charm quarks. In the BBNS ap- 
proach, the charm quarks are treated as hard modes (i.e. rric = 0(m;,)), and therefore they 
are integrated out in the first matching step and do not appear as degrees of freedom in 
SCETi. Alternatively, one may take the point of view that ml ~ Am?,, i.e. rric <C rrih in the 
heavy quark limit. Still, since charm quarks do not appear as external partonic degrees 
of freedom in charmless non-leptonic B decays, they cannot induce endpoint singularities, 
and can, in any case, be treated perturbatively. The effect of the alternative power- 
counting scheme merely amounts to expanding the hard coefficient functions in terms of 
mc/mb, which corresponds to integrating out the charm quarks in the second matching step 
SCETi SCETii. Via the renormalization group running within SCETi one would also 
resum logarithms In rric/ rrib. Numerically, the effect onto factorizable amplitudes should 
be marginal. In the non-factorizable contributions, on the other hand, the crucial effects 
come from the endpoint divergences and chiral enhancement related to light quarks, and 
not from charm quarks, if rric ^ Aqcd- 

In j3] it has been argued that one should also consider the possibility of non-relativistic 
charm and gluon modes in SCETi. This would correspond to operators of the type 

J [CC] NR , 

where the invariant mass of the cc pair happens to be close to 4m^ which corresponds to 
a light-cone momentum fraction for the quark field of the order u = 1 — 4ml /ml — 
0.6. However, since charm modes do not appear as external degrees of freedom, hadronic 
matrix elements of the above operator in SCETi x NRQCD would vanish. One may ask 
the question, where the effect of charm resonances (i.e. non-relativistic cc bound states, 
J/ip, ■ ■ .) would show up in the effective theory framework. As explained above, as charm 
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quarks do not appear as external states, they can be formally integrated out, resulting in 
a quark determinant with charm quarks in the background of soft and coUinear gluon 
fields. The treatment of charm quarks is thus fully inclusive, and the appearance of charm 
resonances resembles the well-known cases of e^e~ — > hadrons or B —>■ Xsi~^i~: When 
integrating the invariant-mass spectrum over a large enough region, the effect of charm 
resonances provides power-corrections to the inclusive rate. 

In our case, the pion distribution amplitude serves as the "detector" with a "sensitivity" 
(pniu = 1 ~ (f l^'b)- Integrating over all momentum fractions m, the effect of charm res- 
onances translates into power-corrections, which should be attributed to the matching 
coefficients of sub-leading operators in SCETj. One can even perform a numerical esti- 
mate of such charm-resonance effects by using the same phenomenological treatment as for 
h s£~^i~ ISOl EI] • We find that the standard pion distribution amplitude is sufficiently 
broad to wash out the effect of exclusive charm resonances. Therefore it is not clear to us in 
what sense the inclusion of non-relativistic modes in SCETj should lead to an enhancement 
of charm penguin contributions in non-leptonic B decays, as has been argued in Also in 
a recent analysis within QCD sum rules sua unnatural enhancement of charm-penguin 
contributions to non-leptonic B decays is not observed, and the perturbative treatment of 
charm quarks seems to be justified. 

2.5.4 5-body operators 

Two examples of 4-quark operators with an additional hard-collinear gluon are 

[U^mjiEALK] - A^ (39) 

[^^cAt^i.M-M ~ A^ (40) 

With similar arguments as above, the first term matches onto a leading-power 6-quark 
operator in SCETn, which corresponds to the factorizable hard-spectator diagrams in the 
QCD factorization approach. The second term gives a l/nib power-suppressed contribution, 
which serves as one of the candidates to compensate the endpoint divergences parametrized 
by Xh- There are a number of further 5-body structures, which we are not going to discuss 
in detail here. 



2.5.5 6-body operators 

The simplest 6-quark operator in SCETi has the form 

The leading-power matching onto SCETn gives a 6-quark operator of the order A^^ which 
corresponds to the l/mf,-suppressed annihilation graphs in the BBNS approach. 

Another possibility is a 4-quark operator with two additional transverse gluons, e.g. 

[^UT^][qsK]AiAt, X'' (42) 
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which amounts to annihilation with additional radiation of hard-collinear gluons. After 
matching onto SCETn, the minimal quark content refers to an 8-quark operator which 
scales as A^^, and is therefore l/ml suppressed with respect to the leading-power contri- 
butions, and should be related to the terms of order X\ in (j7fj) . 

2.5.6 Summary: Factorizable and non-factorizable SCET operators 

We summarize the examples for SCETi and SCETn operators in TableEl The identification 
of the power-corrections as non-factorizable is tentative (without a more detailed analysis, 
which is beyond the scope of this work, we cannot exclude that some of the operators 
identified as non-factorizable at a certain sub-leading power in l/m^ are actually free of 
endpoint divergences in SCETn). The factorization of the two structures contributing at 
leading power (as indicated by ^J), can be understood by applying the same rules as in the 
B iriu case |5] (this point has also been realized in Q). 

Notice that another set of possible operators is obtained by changing one gluon field 
into a photon field. The contributions of such operators are suppressed by the ratio of 
electromagnetic and strong coupling constants, aem/cts- They may be important in isospin- 
breaking observables for B — > ttK decays, if they enter with large Wilson coefficients 
and large CKM elements, and if the l/mb suppression is compensated by large numerical 
factors. 
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Table 2: An incomplete list of the field content of SCETi operators, contributing to 
non-leptonic B decays into light hadrons. See text for details. 
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2.6 Non- factor izable effects from long-distance penguins 




(a) As (b) As As 

Figure 9: Examples for (a) short- and (b) long-distance penguin contributions to a 4-body 
Greens function in SCETj via hard, hard-collinear or soft quark loop. 

It is often argued that a main source of non-factorizable effects m B ^ PP decays 
should be attributed to long-distance penguin topologies [211 122] , where two quarks of one 
of the four-quark operators Oi_6 in the electroweak hamiltonian are contracted to a loop 
which radiates one or more gluons. As explained above, in the effective theory approach, 
long-distance modes in SCETj are hard-collinear and soft quark and gluon fields, i.e. the 
effect of h (and usually also c) quarks is already absorbed into matching coefficients, see 
Fig. Efa). Examples for long-distance penguin diagrams in SCETj are shown in Fig. IHl^b). 
They reflect contributions to one of the 4-body operators (fHU)) discussed above, which 
gives power-suppressed contributions to the decay amplitudes. Remember that above, we 
identified this type of operator as the source for annihilation contributions. Actually, in 
the non-perturbative region the diagrammatic language in terms of soft quark and gluon 
propagators is not appropriate anymore, and the distinction between annihilation and 
penguin contractions is not clear cut. On the other hand, in SCET one has an additional 
correlation between flavour quantum numbers and the type of fields (hard-collinear or soft). 
For the above example, the operators 

fc/it,][4^^hc]/=o and [qshv\[l-^i)ic]i=i , (43) 

where the total isospin of the two hard-collinear fields is different, can be distinguished. 
Notice however that, in general, interactions between soft and hard-collinear quarks can 
lead to mixing of the two operators in SCETi. Notice that by momentum conservation 
both hard-collinear fields are in an endpoint configuration 

= Enl + Ap'^, |Aj9''| ~ A 

= E< + Ap'', |Ap''| ~ A , (44) 

where E ~ mb/2 denotes the pion energy in the B meson rest frame, and n_ and n+ are 
light-cone vectors which satisfy n'^_ = n\ = 0, = 2. The invariant mass of the gluon 

pair is of order y/Arrih, and therefore the light-quark thresholds in the loop give rise to an 
imaginary part in the soft loop. We emphasize that for this configuration the factors of 
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as belonging to the soft gluon and to the endpoint-gluon do not count as perturbative. 
Therefore the latter reflects a mechanism to generate non-perturbative strong phases. On 
the other hand, as already mentioned, this contribution is power-suppressed by at least 
1/ml, in accordance with the general arguments in and and, in principle, the 
endpoint-configuration should be suppressed by Sudakov form factors. On the other hand, 
we could have a sizeable numerical enhancement from chiral factors, but this cannot be 
quantified in a satisfactory way at the moment. 

2.6.1 Scenario 3: Dominance of strong phases in ylc(0, 1/2) 

Sometimes the main source of non-factorizable effects in B ^ tttt decays is attributed to 
matrix elements of the operators ^ [db] [cc] in the electroweak hamiltonian (see e.g. j2Il 
Wl\ for the general argument, and jSSj for a recent hadronic model). This would correspond 
to long-distance penguin topologies with charm quarks in the loop. In the effective theory 
approach this would require to treat 4m^ < Am^. Also, the line-of-reasoning in |lj (see the 
discussion about non-relativistic charm modes, above) has lead the authors to conclude 
that charm penguins are the primary source of strong rescattering phases. 

In a more general setup, such a situation can be simulated by setting the phases 
(/)„(0, 1/2), 0^(2, 3/2) = 180°. Furthermore, we require moderate values for the non- 
factorizable parameters, r„(0,l/2) < 0.5 and r„(2,3/2) < 0.2. Notice that, in general, 
one cannot decide whether the remaining contributions to A^^(0, 1/2) are to be identified 
as long-distance penguin or annihilation topologies (and therefore also the phenomenolog- 
ical discussion in [3£j is covered). The result of this scenario (where the value of the form 
factor is again fixed as 0.26) is shown in Fig. ^1 We observe that the scenario is very 
restrictive, leading to rc(0, 1/2) ~ 0.15, and 0c(O,l/2) ~ 100°. This implies that, based 
on phenomenological assumptions, the theoretical predictivity compared to the more gen- 
eral case in scenario 2 is improved. However, these constraints, used in a CKM analysis, 
probably represent a large theoretical bias. 

We should mention that the phenomenological scenario constructed in j3] differs from 
the one discussed here in two points: (i) In j3] the remaining contributions to ru{I ■, A/) 
are attributed to factorizable hard-scattering contributions, where the relative weight for 
different isospin amplitudes is fixed but the absolute normalization is kept arbitrary, (ii) 
The form factor value is left as a free parameter. In this case the authors found a very 
small value Ff^^"^ ~ 0.17 which seems to contradict the findings from QCD sum rules 
and would point to large hard-scattering corrections relative to the naively factorizing 
terms. We understand this as an indication that the neglect of other non-factorizable 
corrections in ^ is questionable. Since the transition form factor itself and the factorizable 
corrections actually represent one part of the calculation where we have some theoretical 
control on the parametric errors, we believe that one should stick to the present theoretical 
estimates, unless a more direct experimental determination of these terms (from B -kIv 
and B — > 7^//) gives more reliable numbers. 

What remains true is that the assumption about the dominance of strong rescattering 
phases in the amplitude v4c(0, 1/2) cannot be excluded with present data. In view of the 
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already mentioned QCD sum rules result |S2] and the general arguments in Section I2.5.3t 
it does not seem very plausible that these effects are due to charm penguins alone. 
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Figure 10: Comparison of values for a random sample of non-factorizable parame- 
ter combinations, using 0„(0, 1/2) = 0„(2,3/2) = 180°, F^^^ = 0.26, r„(0, 1/2) < 0.5, 
ru(2,3/2) < 0.2 (Scenario 3 - dominance of strong phases from Ac{0, 1/2)). 



2.6.2 Scenario 4: Equal strong phases in A„(0, 1/2) and A^^O, 1/2) 

As explained above, in the standard framework, long-distance penguins would only involve 
light quarks. Imaginary parts are related to the soft momentum regions in these light-quark 
loops. If this was the leading mechanism to generate strong phases, one may indeed expect 
that the corresponding contributions to Au{0, 1/2) and Ac{0, 1/2) differ in moduli, due to 
different Wilson coefficients, and short- distance coefficient functions in SCET, but yield 
the same strong phases 0^(0,1/2) = 0c(O,l/2). (More precisely, this situation is always 
realized if the non-factorizable contributions are dominated by only one operator in the 
AI = 1/2 SCETi hamiltonian, e.g. the operator in Strong phases can be induced 

either by soft quark loops in penguin diagrams, or by soft gluon rescattering in annihilation 
diagrams. As already mentioned, to disentangle the two possibilities does not make sense 
for non-factorizable operators.) The corresponding 4-quark operators with the largest 
Wilson coefficients are 0^2i and therefore it is also conceivable that r„(0, 1/2) > rc(0, 1/2). 
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As one can see from the general set-up in scenario 2 (see Fig. |HI), both these assumptions 
seem to be in hne with experiment. 

We therefore define another scenario 4 with four constraints, and thus only three free 
parameters, ^4^^^(2,3/2) = 0, 0^(0,1/2) = 0„(0, 1/2), and F^^^ = 0.26. The comparison 
with the experimental data is shown in Fig. ^2 Not surprisingly, the largest restriction 
concerns the value of 0„(0, 1/2) = 0c(0, 1/2) which is tuned to values around 150°. But 
the values of r„(0, 1/2) and rc(0, 1/2) which lead to a good description of the data are still 
rather generic, r„(0, 1/2) ~ 0.4 - 0.8 and rc(0, 1/2) < 0.2. 

2.7 Lessons from B ^ tttt 
We conclude: 

• The present data on B ^ Tin decays require non-vanishing non-factorizable correc- 
tions which should be related to l/m^, corrections in the factorization formula (jT]). 
(Here we assume that the CKM elements are given by their SM fit values.) 

• The dynamical origin of these corrections remains a theoretical challenge, and differ- 
ent phenomenological assumptions can accomodate the data. This includes scenarios 
in the spirit of BBNS, where non-factorizable corrections can still be moderate for 
certain hadronic input parameters. On the other hand, the central values of exper- 
imental data seem to point to rather large non-factorizable contributions, in partic- 
ular for the isospin amplitude Au{0,l/2), which can only be reached by pushing the 
hadronic parameters in the BBNS approach to the limits. 

• Different assumptions about the dominance of certain decay topologies are consistent 
with the data. However, the additional assumptions and constraints may lead to a 
strong theoretical bias when used in CKM fits. Whereas it seems safe to neglect 
the non-factorizable contributions to the isospin amplitude ^4^(2, 3/2), in general, 
both large contributions to ^^^(0, 1/2) and Ac{0, 1/2) should be taken into considera- 
tion. They may be related to either long-distance penguin or annihilation topologies 
which correspond to matrix elements of power-suppressed operators in soft-coUinear 
effective theory. 

• The picture of large non-factorizable penguin contributions still appears to be at- 
tractive from the phenomenological point of view. Although they formally appear 
on the level of power-corrections, their numerical impact remains a matter of debate. 
We have given some arguments that in the effective theory framework, one would 
need to consider more complicated diagrams than those taken into account in P or 

in order to cover such effects. We also saw that the distinction between annihi- 
lation and penguin topologies is not clear cut for non-factorizable contributions. We 
also find it unlikely that large non-factorizable effects to ^4^(0, 1/2) are coming from 
charm-quark loops. 



25 



10 



r„(0, 1/2) 


• " '. ' 

. •■' 


. 

.V' ii- ■. 
-■ f' 


. • i'"'' * • •• 




1.1, V-*. ■•. 


i" 






















■ --^-^ 




i" « 




,^\" 





).2 



0.4 0.6 

rc(0, 1/2) 



).2 



10 



50 



100 



0.4 0.6 

0^(0,1/2) 

7*. 
* • 

150 



o.s 



200 



250 



300 



350 



Figure 11: Comparison of values for a random sample of non-factorizable parameter 
combinations, using F^^^ = 0.26, 0„(0, 1/2) = 0^(0, 1/2) and r„(2, 3/2) = (Scenario 4 - 
equal phases for ^^^(0, 1/2), and ^^^(o, 1/2)). 
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• In the extreme case, we may assume that the non-factorizable effects required by 
experimental data are dominated by one non-factorizable operator in the A J = 1/2 
SCETi Lagrangian, leading to equal strong phases for A^^{0, 1/2) and A^^{0, 1/2). 
Comparison with experimental data shows that these long-distance contributions to 
the decay amplitudes, are about an order of magnitude larger than the factorizable 
penguin contributions contained in (0). We note that in the BBNS approach this 
scenario can be simulated by allowing for significantly larger non-factorizable anni- 
hilation contributions X^- 



3 Some Remarks on B Kn Decays 

In the future, we expect new and/or more accurate data within the whole class oi B ^ PP 
decays from the B factories and hadron machines. This is crucial for a better understanding 
of the strong dynamics in charmless non-leptonic B decays. In this work, we limit our 
analysis to some general remarks. We plan to give a comprehensive analysis of -B ^ Kn 
decays in a forthcoming paper. 



3.1 Isospin decomposition for B jtK 

The isospin decomposition of 5 ^ nK amplitudes reads 





HeS 


B-) 




[A„(l/2, 0) + A„(l/2, 1) + A,(3/2, 1)] + 
[A,(l/2, 0) + A,(l/2, 1) + A,(3/2, 1)] , 


(45) 




HeS 


B') 




[-A„(l/2, 0) - A„(l/2, 1) + 2A„(3/2, 1)] + 
[-A,(l/2, 0) - A,(l/2, 1) + 2A,(3/2, 1)] , 


(46) 


{tt+K- 


\Hcs 






[-A„(l/2, 0) + A„(l/2, 1) + A„(3/2, 1)] + 
[-A,(l/2, 0) + A,(l/2, 1) + Ae(3/2, 1)] , 


(47) 




\HeS 




^ Ai^) 


[A„(l/2, 0) - A„(l/2, 1) + 2A„(3/2, 1)] + 
[A,(l/2, 0) - Ae(l/2, 1) + 2A,(3/2, 1)] , 


(48) 



reflecting the isospin relation 

{■K-K^\H^s\B-) + V2{7T°K-\H^s\B-) - {7i+K-\H,s\B°) - V2 {-n^ K^\H^q\B^) = 

For comparison, we also quote the connection with the parametrization used in 

P = A('^)(A,(l/2,0) + A(l/2,l) + A,(3/2,l)) , 



^3/2 



i4>a 



3*0 



(-KM 



(KM 



A^(l/2, 0) + A„(l/2, 1) + A„(3/2, 1) 
A,(l/2,0) + A,(l/2,l) + Ae(3/2,l) 

3A„(3/2,1) 

A,(l/2,0) + A,(l/2,l) + A,(3/2,l) 



. 

(49) 



(50) 
(51) 

(52) 
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eye 



2A„(1/2,1) + 2A„(3/2,1) 
A,(l/2,0) + A,(l/2,l) + Ac(3/2,l) 

1 ^c(3/2,l) 



gee 



eKM A«(3/2, 1) ' 
1 Ae[l/2,1]+A,[3/2,1] 



(53) 
(54) 

e;,M A„[l/2,1]+A,[3/2,1] ^^^^ 

Thus, there are eleven independent isospin parameters for the Kit mode. At the moment, 
there is not enough data available to fix them all independently. 

The factorizable part of the isopin amplitudes in the QCD factorization approach can 
be expressed using the corresponding parameters a^j and Oj^i (the latter coefficients are 
restricted again to the heavy quark limit) as given in pp. 



(3/2,1) 
^r(3/2,l) 
(1/2,1) 

<(1/2,1) 
^r(l/2,0) 

^^(1/2,0) 



A 



ttK 



2 

~Y~ 

AnK 
T2" 

AnK 



K c c \ , -^K-K / , \ 

-r^ Og - a^o ) + (-07 + ag) 



Aktt 



2ai - 3r^ < - Sa^^Q) + ^ (-2a2 + 3a7 - 309) 



6 



K c c \ I ^Ktt / N 

-r ag - flio + (07 - ag) 



-2ai - 3r^ < - 3<o) + 



A 



K-K 



-2a2 — 3a7 + 309) 



2ai + 4a^ + <o + 4r J a;^ + 



(56) 
(57) 
(58) 
(59) 
(60) 
(61) 



3.2 SU{?))f connection between B ^ tttt and B ttK 

A well-defined procedure to reduce the number of independent hadronic amplitudes in 
the analysis of non-leptonic B decays is to use the limit of SU{3)f fiavour symmetry. 
For a comprehensive summary see Appendix A of jSZl-'' Together with the structure of 
the effective hamiltonian, SU{3)f symmetry results in the following relation between the 
isospin amplitudes for B ^ tttt and B —>■ ttK decays 



Ap(2,3/2) = A,(3/2,l) 



P 



u, c 



(62) 



These equations represent in principle four additional relations between the real parameters 
entering the definition of the isospin amplitudes, namely two phase relations and two ones 
connecting the moduli. Using the parametrization ()10I12|) and the approximate relation 
(jni) for B ^ TTTT amplitudes together with the parametrization for B Kir amplitudes 
quoted above, Eqs. translate into 

3 Cg + Cio 



q e 



EW 



2eKM Ci + C2 



0.69 



(63) 



''Note that we use another sign convention for the definition of isospin amphtudes 
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and 



1^3/2 P| _ l(l + xe-^)T| 
u(^)| ~ uW| • 

I All I 1^" I 

The first equality states that in the 5*^7(3) p limit the parameter q e'^^ is a pure short- distance 
quantity, i.e. a real number that can be calculated in terms of Wilson coefficients and CKM 
factors. This is a well-known result where only the structure of the electroweak effective 
hamiltonian enters ^Tj QHl- The second relation connects the overall normalization, i.e. 
the dominant amplitudes T for B ^ titi, and P ioi B ^ ttK in the naive factorization 
approach. These are known to receive sizeable SU{3)f breaking corrections, already in 
naive factorization through the ratio 



[m: 



2^ 



1.2-1.3. (65) 



There is no further exact SU (3) p connection between B ^ titi and B — > Ktt decays, 
because the SU (3) p multiplets involve also other decays, namely B — * KK and Bs Kn 
decays for AS" = and Bg KK and B^ mi decays for AS* = 1. However, using 
additional information by anticipating some experimental data, one can derive further 
SU{3)f relations between the 5 — >■ tttt and B — > Kir modes. 

If one assumes that A[B^ K^K^] (which in the SU{3)f limit corresponds to 
A[B^ TT+TT^] —>■ 0) SU{3)p symmetry implies the following additional relations between 
isospin amplitudes 

Ap{0, 1/2) = Apil/2, 0) - Ap(l/2, 1) , p = u,c (66) 

In the factorization approach this corresponds to neglecting exchange/annihilation topolo- 
gies. It can be easily translated into four additional relations conecting two phases and 
two moduli in the vrvr and in the Kir mode. 

Often some further empirical assumptions are used within the SU{3)f analysis, namely 
Adir[B^ rc^K'^] = or BR[B^ K+K'^] = 0. This allows to eliminate two more 
parameters within the SU{3)f analysis, namely (see (jSHI)) 

ea e''^^ = . (67) 

Clearly, a complete analysis of the systematic error of such a procedure has to take into 
account an estimate of possible SU{3) f breaking effects, and of the uncertainties related to 
the phenomenological/empirical assumptions. The simplest approach to estimate flavour- 
symmetry violating effects, is to restrict oneselves to factorizable amplitudes as in 
However, in case of large non-factorizable QCD effects, which appear to be necessary 
to explain the B —>■ tttc data, one consequently has to take into account the possibility of 
sizeable flavour-sjTumetry violating effects for such contributions, too. Naively factorizable 
SU{3) F breaking effects thus can only provide an order-of-magnitude estimate of flavour- 
symmetry breaking in general. 
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In the BBNS approach, apart from different decay constants and form factors for pions 
and kaons, SU{3)f breaking enters through the different (moments of) hght-cone distri- 
bution amphtudes. As we aheady pointed out in Section 2, the universal treatment of 
non-factorizable parameters Xh and Xa in all B — > PP isospin amphtudes is question- 
able, because, on physical grounds, we expect the related low-energy dynamics to depend 
on the light quark masses in an essential way. In this context we note that already in 
the framework of the BBNS approach, there are additional sources of SU{?))f breaking. 
According to jHH]! one has to consider strange-quark mass corrections to the twist-3 am- 
plitude 0p (which in the BBNS approach parametrizes the endpoint behaviour of the 
non-factorizable diagrams). We have (neglecting again all other terms that vanish in the 
asymptotic limit) 

= 1 + pI,k {-\ Cl'\2u - 1) - ^ Cl'\2u - 1)) , (68) 
where C^(^) are Gegenbauer polynomials, and 

2 

p\ = ^ = C(m,/A) ~ 5 - 10% ^ . (69) 



K 



In particular, at the endpoint, m — 1 we obtain 

1 - ij^-l ^ 15 - 35% (70) 

which is the typical size of flavor symmetry breaking usually assumed in other phenomeno- 
logical analyses. Notice that the above considerations point to smaller non-factorizable 
effects for kaons than for pions, which is different from the behavior of the factorizable 
terms, where fx > f-w and F^^^ > F^^'^ . As this estimate does not include any dynam- 
ics related to the actual strong rescattering, we cannot exclude that the true effect might 
be even larger. 

Finally, we note the parameter qcC^^" in is neglected in many 5'[/(3)f analyses (see 
for example JI]), using the fact that in the limit of naive factorization, this parameter is 
colour-suppressed. However, also the parameter x within the i? ^ tttt is small in naive 
factorization due to this argument, but nevertheless that parameter is found to be large 
when compared to experimental data (see section 2). 

Thus, we conclude that in view of large strong phases and large contributions from 
colour-suppressed topologies observed in experiment, any other prediction, following from 
the factorization approximation, should be critically re-analysed. 



3.3 i?„-puzzle 

The present data on the Kt^ branching ratios can be expressed in terms of three ratios 
P _ BR[i?° ^ n-K+] + BR[^° ^ n+K'] _ , 
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1 BR[i?° ■ 

2 BR[Bj 



Tr-K+] + BR[Bl 



TT+K-] 




(72) 



Rc 



BR[i?+ 



1.16 



+0.13 
-0.11 



(73) 



This result appears somewhat anomalous, when compared, for example, with the approxi- 
mate sum rule proposed in [1211111113 which leads to the prediction Rc = Rn, or from the 
comparison with the vrvr data using the SU{3)f symmetry approach. In particular, the 
ratio Rn appears to be smaller by about two standard deviations than one would expect. 
It is important to note that the deviation of i?„ and Rc from unity is solely due to isospin- 
breaking effects |^ . The amount of short-distance isospin breaking in the Standard Model 
is too small to explain the experimental number. Whereas the authors of ^1] argue that 
this may point to an interesting avenue towards new physics in electroweak penguin oper- 
ators, the collaboration in |12j considers this deviation as a statistical fluctuation, which 
is consistent with the Standard Model - even when the generalized SU (3) p constraints, 
and (jnni), are enforced. Not surprisingly, the B Kit data has triggered several 
new-physics analyses (for the very recent literature see [TTl liH} li^ E^l 1^). 

Apart from the theoretical questions about the interpretation of the data, there are 
also two experimental issues which have to be clarified: 

• Radiative corrections to decays with charged particles in the final states may not 
have been taken into account properly in the experimental analysis, an effect which 
is expected to lead to an increased branching ratio of these modes ^2] and which 
could bring Rn closer to unity. 

• It has also been argued in [HU] that the present pattern could result from underesti- 
mating the 7r° detection efficiency which implies an overestimate for any branching 
ratio involving a 7r°. The authors of [50j propose therefore to consider the ratio 
{RnRcY^"^ in which the 7r° detection efficiency cancels out. Again, the experimental 
value for this quantity is closer to unity. 

But even when the experimental situation is clarified and the experimental accuracy 
will be significantly improved, one has to critically analyse whether the data on B ^ mi 
and B nK decays point to new physics (for example to isospin breaking via new degrees 
of freedom as discussed in jSHES); or whether they can be explained by non-factorizable 
SU (3) F- or isospin-violating QCD and QED effects within the SM. Obviously, the inclusion 
of more hadronic parameters to include such effects can only improve the phenomenological 
situation, compared to the analyses based on flavour-symmetry or the constrained BBNS 
scenario. On the other hand, the order-of-magnitude of flavour-symmetry breaking effects 
should be consistent with the theoretical and experimental estimates of the non-factorizable 
contributions in the B — »• tttt sector. 

In order to illustrate this point, we discuss a toy model which is inspired by the result of 
the previous section. We start with the assumption that all nonfactorizable contributions 
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are saturated by long-distance QCD and QED penguin contributions. To estimate the pos- 
sible numerical effect, we multiply the factorizable QCD and QED penguin by a common 
enhancement factor (1 — APnf), where APnf is a complex parameter. Within the BBNS 
approach (see (fT^ and (jSEI)) the QCD penguin function pQ'^^ appears in the combination 
+ ag r^, while the penguin function occurs in the combination a^g + of '^x- The two 
SU{3)p constraints, (jU^ and (jUBj) . are both compatible with this procedure by construc- 
tion. From the B ^ titt decays, we find that the long-distance QCD penguins require an 
enhancement by about an order of magnitude (long-distance QED penguins are negligible 
in S — >■ tttt). Consequently, the long-distance QED effects in 5 — > irK are enhanced by 
the same amount, which leads to significant changes in the isosp in- violating parameters 
ge*"^ and gc-e*"^^. We find typical values in the range 0.3 < g < 0.8, —30° < uj < 30°, 
0.05 <qc < 0.3, and -150° < coc < 150°. 

More generally, in our framework the expansion parameter that suppresses long-distance 
isospin-violating effects in B ^ ttK is given by 

r,(0, 1/2) . ~ 5 - 10% , 

where the latter number refers to typical numbers for rc(0, 1/2) = 0.1 — 0.2 observed in 
B 7171. This seems to be in the right ball park to at least partly explain the deviation 
of Rn from unity. 

If realized in nature, the latter scenario should also have some impact for other "puz- 
zling" B decays: In B ^ 4'Ks the apparent hierarchy r„(0, 1/2) > rc(0, 1/2) in B 
decays may translate into a moderate enhancement of the CKM-suppressed penguin am- 
plitude. This may lead to some deviation of the extracted value of sin 2j3es from the value 
found in B ^ J/ijjKs, but is probably not sufficient to explain the present BELLE mea- 
surement jSH] for this quantity. Non-factorizable penguin contributions can be even more 
important in B ^ rj'K decays, because for decays into singlet mesons some non-factorizable 
operators appear at lower order in the l/m;, expansion than for the octet ones (see the 
discussion in Section [2.5^ . According to [23 this mainly implies that the theoretical un- 
certainties in such decays presently are too large to extract information on short-distance 
quantities in a reliable way. 

Our main conclusion is a conservative one, namely that the unsatisfactory theoretical 
understanding of non-factorizable power-suppressed effects in charmless non-leptonic B 
decays prevent us from identifying new-physics effects in these observables, at least on 
the level of present experimental significance of certain "puzzles". On the other hand, 
the comparison of different possible approximation schemes, used to reduce the number 
of unknown hadronic parameters, gives a handle to estimate model-dependent systematic 
effects in CKM studies. It may also shed some light on the dynamical origin of non- 
factorizable effects, which may stimulate further studies with non-perturbative methods. 
In particular, a systematic classification of power-suppressed matrix elements from non- 
factorizable SCET operators should give an alternative scheme compared to the traditional 
classification in terms of flavour topologies. 
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